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1  Summary 

This  report  svraimarizes  work  carried  out  under  ONR  Grant  Number  N00014-97-1-0053  at  the  University 
of  Michigan,  and  coordinated  efforts  under  ONR  Contract  Number  N0014-00-C-0057  at  Veridian 
Systems  Division,  Inc.  The  goal  of  this  effort  was  to  improve  the  ability  to  predict  die  free-surface  flow 
in  the  near  field  of  surface  ships.  The  spo^ific  objective  were:  to  develop  modeling  approaches 
appropriate  to  the  near -surface  region  h  turbulent  free  surface  flows  and  to  obtain  data  for  validation  of 
CFD  predictions  of  free-«urface  turbulent  flows.  The  context  of  this  modeling  effort  was  the  R^olds- 
averagedNavier-Stokes  (RANS)  equations,  as  implemented  in  the  code  CFDSHIP-IOWA.  The  specific 
focus  was  the  modeling  of  unsteady,  unresolved  mjrface  waves  ^nerated  by  the  near-field,  turbulent  flow 
around  a  surface  ship ,  including  die  generation  and  propagation  of  these  waves.  The  main  results  of  fliis 
study  can  be  summarized  as  follows:  1)  The  fiamewoik  outlined  in  Hong  &  Walker  (2000)  was  shown  to 
be  reasonable  for  the  analysis  and  prediction  of  turbulent  firee-surface  flows.  2)  The  surface  current  is 
caused  by  a  combination  of  turbulmt  surfece  fluctuations  and  turbulence  anisotropy,  widi  surf^e 
fluctuations  dominating  in  hi^Froude-number  flows,  and  anisotropy  in  low-Froude-number  flows.  3) 
The  surface  fluctuations  effecte  can  be  approximated  by  a  modified  stress  boxmdary  condition  where  the 
shear-stress  is  related  to  the  gradient  of  the  surface-fluctuation  variance.  4)  A  computational  approach  for 
predicting  the  production  and  evolution  of  turbulence-generated  waves  was  developed  and  the  resulte 
compared  favorably  to  experimental  data.  5)  Surface  fluctuation  measurements  were  carried  out  for  a 
firee-surface  jet  flow  and  probability  density  function  for  tiie  surfece  elevation  was  shown  to  be  relatively 
Gaussian.  ^  A  method  for  measurement  ofthe  directional  wave  spectrum  was  developed  and 
implanented  and  the  directional  spectnan  of  turbulaxce  generated  waves  was  characterized  for  a  fra- 
surfacejet. 


2  Background 

The  results  of  Walker  et  al.  (1995)  for  turbulent  free-surfece  jets  show  that  there  is  significantly  more 
wave  generation  in  high-Froude-number  flows  than  in  low-Froude-number  flows.  There  is  also  a  marked 
reduction  in  the  turbufence  kinetic  energy  levels  in  flows  at  hi^  Froude  number;  the  implication  being 
that  the  turbulence  kinetic  aier^  is  being  transferred  to  wave  energy,  which  then  propaptes  away  in  the 
form  of  surface  waves.  This  loss  of  energy  (roughly  20%  for  Fr  =  8,  relative  to  Fr  =  1)  will  affect  the 
ultimate  downstream  evolution  of  the  flow.  Observations  ofthe  near-field  flow  around  both  fulhscale 
and  model-scale  ships  reveal  si^ficant  turbulent  wave  geno’ation  in  these  flows.  As  a  result,  accurate 
prediction  of  these  flows  requires  that  the  generation  of  waves  by  turbulence,  as  well  as  their  subsequent 
propagation,  be  modeled  appropriately.  Since  these  waves  are  not  resolved  either  spatially  or  temporally 
in  Reynolds -averaged  calculations,  their  effects  also  need  to  be  included  in  the  turbulence  modeling  for 
the  sub-surface  flow. 

RANS  codes  such  as  CFDSHIP-IOWA  (Tahara  &  Stem  1996)  can  model  a  fully  deformable  mean  free 
surface,  and  thty  impose  averaged  versions  of  the  kinematic  and  dynamic  fr^-^ir&e  boundary 
conditions  at  fee  mean  free  surface.  Careful  examination  of  fee  instentaneous  equations  and  boundary 
conditions,  however,  indicates  feat  this  approach  is  incompatible  wife  unsteady  surface  fluctuations,  and 
is  ftmdamentally  incorrect  when  such  unresolved  surface  fluctuations  are  present.  This  is  because  fee 
boundary  conditions  apply  at  fee  instantaneous  position  of  fee  free  surface,  and  averaging  fee  equations 
does  not  allow  fee  conditions  to  be  applied  at  fee  mean  free  surface  inst^d,  except  in  fee  case  where  the 
mean  and  instantaneous  free-^urface  locations  are  the  same  (i.e.  no  fluctuations). 
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Including  the  unresolved  surface  fluctuations  in  the  Reynolds  averaging  process  leads  to  a  set  of 
equations  which  ^plies  across  the  mtire  air-\^er  fluid  volume  (Hong  &  Walker  2000).  One  form  for 
these  equations  is 


dt 


(1) 


where  (t/,) 


dxj 


du/ 

dx,i 

-p{uiUj)+pgi-tf(^z) 
K 


= pU  J p  is  a  density-wei^ted  (Favr6)  average,  and  the  overbar  indicates  a  time 


(2) 


average;  p. 


pand  /are  density,  viscosity  and  surface  tension,  respectively,  and  ^  is  an  ava-age  measure  of  the 
surface  curvature.  Here,_/(z)  is  the  p.d.f.  of  the  firee-surface  fluctuations,  Po  and  Po  are  the  water  daisity 
and  viscosity,  and,  for  a  single -valued  free  surface. 


P W = -  ll/W  dz^  p(2) = p^^l  -  J^/(2)  ife j 


(3) 


at  a  given  position.  The  equation  govaning  the  mean  free^urface  position  is 


dt 


+ 


M  =  p 


which  is  derived  from  the  time-averaged  form  of  the  continuity  equation.  Here  G  =  0  defines  the  free 
surface.  It  should  be  noted  that  the  set  of  equations  (1),  (2)  and  (4)  reduce  to  the  levetset  method  of 
Chang,  et  al.  (1996)  in  the  limit  of  laminar  flow. 


Solution  of  the  above  eqxiations  suffers  from  the  same  closure  problems  as  the  conventional  RANS 
equations,  with  the  additional  unknown  of  the  p.d.f.  of  the  surface  elevation  fluctuations.  In  addition  to 
the  viscous  terms  associated  with  the  mean  strain  rate,  there  is  also  an  additional  viscous  term  caused  by 
the  fluctuating  strain-rate  gradients;  both  of  these  are  expected  to  be  small  in  high-R^olds-number 
free-surfece  flows.  Assumption  of  a  suitable  functional  form  for  j^z)  (e.g.  Gaussian)  reduces  the  closure 
problem  to  one  of  specifying  flie  m^  and  variance  of  the  surfece  elevation  distribution  (in  addition  to 
the  specification  of  the  Reynolds  stresses).  The  mean  surface  elevation  is  known  from  (4),  file  surface 
elevation  variance  requires  additional  information. 

If  the  surface-elevation  fluctuations  are  driven  primarily  by  the  local  sub-surface  pressure  field,  then  the 
local  surface-elevation  variance  can  be  modeled  in  terms  of  the  local  turbulence  properties.  If,  however, 
the  main  contribution  to  flie  surface  fluctuations  is  in  the  form  of  waves,  then  the  propa^on  and 
refraction  of  those  waves  must  be  modeled  accurately  in  order  to  predict  the  surface  fluctuations.  The 
oceanographic  community  h^,  in  recent  years,  developed  statistical  approaches  to  modeling  flie 
propagation  of  surface  waves.  These  thirdrgeneration  wave  models  (Komen  etal.  1996)  describe  the 
evolution  of  the  wave  spectrum  in  spm:e  and  time,  and  in  so  doing,  retein  the  wave-like  nature  of  the 
surface  disturbances.  They  account  for  advection  of  wave  enargy  by  the  fluid  velocity  as  well  as  the 
wave  group  velocity  (radiation  effecte),  and  interaction  of  the  waves  witii  velocity  gadients  (refiaction, 
wave-c\ment  interactions).  Since  wave  energy  is  not  a  ‘scalar’  quantity— each  bit  of  energy  has  a 
proj^^tion  qjeed  and  direction  associated  with  it — ^this  fype  of  wave  spectrum  model  is  the  lowest-order 
statistical  model  which  can  be  ^plied  with  the  ejqiectation  of  obtaining  realistic  results  for  turbulence¬ 
generated  waves. 
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The  action  balance  model  is  expressed  as 


dUj  dN 

dXi  dkj 


-vk^N+SiXt,ki); 


(5) 


here,  N  -  N(ki^i)is  the  wave  action  spectrum  which  varies  in  space  and  N  =  E/a(k),  where  E  is  the  energy 
spectral  density  of  the  wav^  and  o  is  the  intrinsic  frequency  of  the  waves  for  a  given  wavenumber,  as 
determined  from  linear  wave  theory.  Q  is  the  vector  group  velocity  of  the  waves,  and  5  is  a  source  term 
representing  turbulent  wave  generation,  or  other  effects.  The  second-from-last  term  in  (5)  represents 
viscous  dissipation  of  wave  energy,  which  impacte  primarily  the  short  waves. 


3  Results 


The  objective  of  this  study  was  to  develop  an  approach  to  predicting  turbulent  free-surface  flows  in  the 
context  of  the  conventional  Reynold-averaged  Navier— Stokes  equations.  In  what  is  presented  below, 
the  main  resulte  are  described.  These  include  results  on  physical  understanding  in  free-surface  turbulent 
flo\^,  model  developmait,  and  single -point  surface-elev^on  measuremente,  as  well  as  wave  spectrum 
measurements.  These  measurements  can  be  used  to  ^de  development  of,  and  serve  to  validate, 
modeling  approaches  for  turbulence-generated  waves. 


The  surface  current,  a  large  lateral  which  develops  velocity  n^  the  free  surface,  is  a  mqor  feature  of 
many  inhomogeneous  turbulent  free-surface  flows  (see  Walker  1997  for  a  review).  The  new  Reynolds 
avoided  equations  (2),  above,  forms  the  basis  for  a  new  analysis  vriiich  i^allels  tiiat  of  Walker  (1997), 
but  includes  tire  effects  of  unsteady  surface  deformations.  This  more  general  analysis  results  in 


-  3F  ■'dV  -  3F 
ax  ay  dz 


4z‘ 


dn'" 

dy  ’ 


(8) 


where  rf  is  the  r.m.s.  surface  fluctuation  level.  This  indicates  that  die  surface  current  is  driven  by  the 
lateral  gMient  of  either  the  Reynold-stress  anisotropy  or  the  surface  fluctuation  variance.  In  high 
Froude  number  flows  where  the  anisotropy  is  small,  the  surface  fluctuation  term  dominates,  while  in  low- 
Froude-number  flows,  the  surface  fluctuation  variance  is  small,  but  the  anisotropy  is  large.  The  details  of 
the  anatysis  leading  to  the  above  result,  along  a  demonstration  that  the  result  is  consistent  with  available 
experimental  jet  data  appears  in  Hong  &  Walker  (2000). 


Surface  fluctuation  efftets  can  only  be  included  in  non-level-set  RANS  codes  in  an  approximate  fa^on. 
For  small  surface  fluctuation  levels,  the  surface  fluctimtion  terms  in  (2)  can  be  d^oupled  from  the  field 
equations  and  an  effective  dynamic  free-surface  boundary  condition  can  be  derived  which  include  the 
effects  of  surface  fluctuations.  The  resulting  r^proximate  boundary  condition  is 


Txz  =4 


^dU  dW 
ydz  ^  etc 


1  drt 


A 


dx 


%  - 


iL  iE 

dz^  dy 


—  2  PS 


.'2 


dy 


Here,  the  sur&e  fluctuations  are  represented  as  an  applied  shear  stress  acting  on  the  m^  free  surfece. 
Details  of  the  derivation  are  provided  in  Walker  (2000), 


A  computational  approach  was  developed  in  which  the  exact  equations  derived  by  Hong  &Walker  (2000) 
are  reduced  to  an  qjproximate  form  for  small  surfece  fluctuations.  This  ^proach  was  described  in 
Walker  (20006, 2001a)  and  allows  hi^-Froude-number  flows  with  imsteady  surface  fluctuations  to  be 
treated  with  a  conventional  RANS  code  such  as  CFDSHIP  -IOWA  (Tahara  &  Stem  1996),  ratha-  than 
using  a  level-set  ^proach.  Surface  fluctuations  are  modeled  via  coupled  solution  of  a  wave-action 
spectrum  model  their  effect  on  the  mean  flow  are  captured  using  an  ^proximate  stress  condition  at  the 
free  surface.  To  accomplish  this  a  wave-action  balance  solver  was  added  to  the  CFDSHIP -IOWA  code. 
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The  effects  of  the  free  surface  on  the  sub-surface  turbulence  at  low  Froude  number  are  c^tured  using  a 

near-airface  anisotropy  model. 

The  approach  developed  were  applied  to  turbulent  foe-surfece  jet  and  comprisMis  of  the  results  to 

experimental  data  were  made  in  Walker  (2001<i,  20016).  First  a  low-Froude-number  jet  was  calculated, 
where  there  was  little  wave  generation  and  near-surface  anisotropy  daninat^  the  behavior  of  the  flow; 
here  the  ^proach  performed  quite  well  when  compared  to  low-Froude-number  jet  data.  A  one-way- 
coupled  prediction  of  ahigh-Froude-number  j^  was  then  presented,  where  the  surface  fluctuations  were 
preicted  using  the  wave-action  model,  but  die  surface  fluctuations  were  not  used  in  the  boundary 
conditions  for  flie  subsurface  flow  (i.e.  the  subsurface  flow  is  treated  as  a  low-Froude-numbo-  flow).  The 
predicted  surface-fluctuation  distribution  was  in  good  agreement  with  the  observations.  Finally,  a  two- 
way  coupled  calculation  was  presented,  where  the  surface  fluctuations  were  used  in  die  boundary 
conditions  on  the  subsurface  flow.  This  agrees  less  well  widx  the  experiments;  the  magnitude  of  the 
outward  velocity  at  the  free  surface  was  substantially  under-estimated.  This  was  believed  to  be  due 
mainly  to  a  lack  of  a  complete  model,  enable  of  switching  between  the  approximate  stress  condition  in 
regions  of  large  surface  fluctuation,  and  an  anisotropic  near-surface  turbulence  model  in  regions  with 
small  surface  fluctuations. 

Surface-elevation  measuremente  have  been  undertaken  to  quantify  the  fluctuation  level  for  die  free 
surface  above  a  turbulent  jet.  These  measurements  were  accomplished  with  a  sin^e  -point  laser-induced 
fluorescence  technique  where  the  water  is  dyed  with  a  fluorescent  dye,  and  the  water  surface  is 
illuminated  from  above  with  an  Argon-ion  laser,  and  are  described  in  Hong  (2000).  The  fluorescent 
emission  along  the  beam  path  is  observed  with  a  line-scan  camera,  and  the  bright-to-dark  transition 
indicates  the  position  of  the  free  surface.  These  measurements  were  used  to  characterize  the  probability 
daisity  frinction  for  the  free-surface  elevation  as  a  function  of  position  on  the  water  surface.  Comparison 
of  plots  of  the  experimentally  determined  p.d.f.s  were  made  to  Gau^ian  distributions  calculated  using 
the  experimentally  determined  mean  md  standard  deviation  to  establish  that  the  free  surface  fluctuations 
are  nearly  Gaussian.  At  some  locations,  (e.g.  on  the  centerline  far  downstream),  the  mode  of  the 
distribution  is  located  slightly  on  the  positive  side  of  the  mean.  These  results  demonstrate  and  reinforce 
the  fact  that  for  this  flow,  and  others  of  practical  interest,  the  free  surface  position  is  not  well  defined  and 
must  be  treated  statistically. 

The  statistical  description  of  the  turbulence-generated  wave  field  is  the  directional  wave  spectrum, 
iiidic^g  the  energy  density  for  waves  of  different  wavenumbers  (wavelengths)  and  propagation 
directions.  The  evolution  of  the  wave  spectrum  is  governed  by  a  dynamical  equation,  in  which  the 
generation  of  waves  by  subsurface  turbulence  is  represented  by  a  source  term,  and  the  generated  waves 
interact  with  mean  velocity  field,  resulting  in  wave  refraction  and  wave-energy  changes  locally.  A 
mefliod  was  developed  which  allows  measurement  of  the  directional  wave  spectrum.  This  approach  has 
yielded  the  fiist-ever  measurements  of  the  wave  spectrum  of  turbxdence-generated  waves. 

To  determine  the  wave  spectrum,  surface-elevation  me^uremente  over  an  area  sufficient  to  determine  the 
wavelength  of  tiie  waves  present  are  required.  It  is  also  necessary  to  determine  the  propaption  direction 
of  tile  wav^.  A  snapshot  of  the  wave  field  yields  an  estimate  of  wave  direction,  but  there  is  a  1 80  degree 
directional  ambiguity.  To  resolve  this  ambiguity  (i.e.  to  determine  whether  the  waves  are  moving  to  the 
left  or  right),  temporal  information  is  required. 

Directional  wave-energy  spectrum  measurements  were  accomplished  using  the  method  described  in 
Lyzenga  et  al.  (2003,  in  preparation).  This  approach  uses  a  vertical  laser  sheet  to  illuminate  the  water 
sur^e  along  a  line.  The  water  contains  a  fluorescent  fye,  and  the  resulting  fluor^cent  emission,  which 
indicates  the  water  surface  elevation  profile  along  the  line  of  tiie  laser  sheet,  is  captured  using  a  CCD 
camera.  The  laser  is  pulsed  and  paim  of  images,  separated  by  a  diort  time  interval  (  Af « 1  ms  ),  are 
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captured.  A  number  of  image  pairs  are  c^tured  and  each  image  is  then  interrogated  to  determine  the 
instantaneous  one -dimensional  surface  elevation  profile  along  the  line  defined  by  the  laser  sheet.  From 
these  profiles,  the  spatial  cross-covariance  fimction  along  the  line  for  time  lag  At  is  estimated.  The  sheet 
is  dien  rotated  a  few  degrees  about  its  vertical  axis  and  the  procedure  repeated,  building  up  a  polar 
representation  of  file  two-dimensional  spatialcross-covariance  fimction  for  time  lag  At .  The  directional 
wave  spectrum  for  a  given  spatial  location  can  be  calculated  fi^om  this  two-dimensional  cross-covariance 

fimction.  Wave  spectrum  measurements  were  carried  out  for  a  free  suif^e  jet  (Walker  eifll.2(K)2).  The 

me^ured  wave  spectrum  above  the  jet  axis  indicates  that  the  dominant  waves  are  propa^ting  upstream 
agamst  the  jet  velocity,  this  is  consists!  with  the  idea  that  energy  would  be  concentrated  in  waves  tbat 
are  ‘blocked’  by  the  opposing  velocity.  Further  fi-om  the  jet  axis,  it  is  clear  that  die  dominant  wav^  are 
propagating  outward,  away  from  the  jet  axis,  consistent  with  waves  observed  in  the  shadowgraph  image. 
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